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ABSTRACT 

We search for and find fossil evidence that the distribution of the spin axes of galaxies 
in cosmic web filaments relative to their host filaments are not randomly distributed. 
This would indicate that the action of large scale tidal torques effected the alignments 
of galaxies located in cosmic filaments. 

To this end, we constructed a catalogue of clean filaments containing edge-on 
galaxies. We started by applying the Multiscale Morphology Filter (MMF) technique 
to the galaxies in a redshift-distortion corrected version of the Sloan Digital Sky Survey 
DR5. From that sample we extracted those 426 filaments that contained edge-on 
galaxies {b/a < 0.2). These filaments were then visually classified relative to a variety 
of quality criteria. These selected filaments contained 69 edge-on galaxies. Statistical 
analysis using "feature measures" indicates that the distribution of orientations of 
these edge-on galaxies relative to their parent filament deviate significantly from what 
would be expected on the basis of a random distribution of orientations. Fewer than 
1% of orientation histograms generated from simulated random distributions show the 
same features as observed in the data histogram. 

The interpretation of this result may not be immediately apparent, but it is easy 
to identify a population of 14 objects whose spin axes are aligned perpendicular to 
the spine of the parent filament (cos 6* < 0.2). The candidate objects are found in 
relatively less dense filaments. This might be expected since galaxies in such locations 
suffer less interaction with surrounding galaxies, and consequently better preserve their 
tidally induced orientations relative to the parent filament. These objects are also less 
intrinsically bright and smaller than their counterparts elsewhere in the filaments. 

The technique of searching for fossil evidence of alignment yields relatively few 
candidate objects, but it does not suffer from the dilution effects inherent in correlation 
analysis of large samples. The candidate objects could be the subjects of a program 
of observations aimed at understanding in what way they might differ from their 
non-aligned counterparts. 

Key words: Cosmology: theory - large-scale structure of Universe - observations - 
galaxies: formation - evolution - kinematics and dynamics - methods: data analysis 



1 INTRODUCTION 

We are searching for observational evidence that galaxies 
are aligned with the large scale structure in which they are 
embedded. Searches for large scale galaxy alignments have 
a long and chequered history going back many decades. It is 
only with the advent of the great galaxy redshift catalogues, 
2dFRS and SDSS, of accurate systematic galaxy photome- 
try and of effective techniques for segmentation of the galaxy 
distribution into voids, filaments and clusters that we can 
now co nfidently addre ss this phenomenon. The recent re- 
view by ISchafen l|2008l ) covers many aspects of the subject 



of galaxy alignments from both the theoretical and observa- 
tional points of view. 

In this paper we specifically search for direct evidence 
for the alignment of the angular momentum of galaxies in 
the SDSS (DR5) catalogue of galaxies with the filamentary 
features in which they are embedded. 

Our approach in this paper is novel and so we briefiy 
review what has been achieved using other, more standard, 
methods. We also review the theoretical motivation driving 
our approach: it is our assertion that much of the large scale 
systemic alignment (as opposed to mere pairwise correla- 
tions) is driven by the tidal fields of the large scale structure. 
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Much of that motivation is based on large scale numerical 
simulations. 



1.1 The Cosmic Web 



Large galaxy surveys such as the 2 dF (ICoUess et aL 1200 ll ) 
and the Sloan Digital Sky Survey (|York et al. 200C) have 
unambiguously revealed an intricate network of galaxies: 
the cosmic web. The cosmic web can be described as a 
mixture of three basic morphologies: blob-like dense and 
compact clusters, elongated filaments of galaxies and large 
planar walls which delineate vast empty regions referred 
to as voids JZeldovich 1970; Shandarin & Zcldovich 198Si). 
iBond et al.l ( 19961 ) emphasized that this web-like pattern is 
shaped by the large scale tidal force fields whose source 
is the inhomogeneous m a tter d istribution itself (see also 
Ivan de Weveaert fc Bondl (|2009l )). Within the large scale 
matter distribution, filaments are characterized by a strong 
and coherent tidal force field: they are tidal bridges in be- 
tween massive cluster nodes. 



1.2 Galaxy Alignments in the Cosmic Web 

With the large scale tidal force field responsible for the mor- 
phology and structure of the Cosmic Web, we expect an in- 
timate link between the Cosmic Web and other prominent 
tidal manifestations of the structure formation process. Ac- 
cording to the tidal torque theory tidal forces generate the 
angular momentum of collapsing halos. Thus we would ex- 
pect the angular momentum of cosmic haloes to be corre- 
lated with the features of the cosmic web in which they are 
embedded. 

Such correlations would arise from the local tidal in- 
fluence of those structures, though they would also be de- 
stroyed by loca l dynamical processes such as merging and 
orbital mixing l|Couttj|l996h . This suggests filaments to be 
the best places to look for alignments: not only are they sites 
with a strong and coherent tidal force field, but galaxies in 
filaments are also likely to have a less disturbed history than 
galajdes in cluster environments. 



1.2.2 The Sloan Digital Sky Survey and 2MASS 

The Sloan Digital Sky Survey (SDSS) has opened up many 
possibilities for analysing the distribution of galaxy ori- 
entations on a variety of scales. Local, small scale (< 
0.5/i~^Mpc), co rrelations among galaxy orientations have 
beei i reported bylSlosar fc Whitd (l2009l).ISlosar et al.l(l2009h 
and I Jimenez et al] (|2010| ). while iBrainerd et al. I 1I2OO9F and 
ICervantes et al ( 2010l ) performed an extensive correlation 
analysis of the relative alignments of pairs of galaxies in the 
DR7 release of the SDSS. The latter authors detected a sig- 
nificant correlation out to projected distance of 10 Mpc 
and were even able to show differences with gala^cy bright- 
ness. The result compared well with the same analysis on 
simulated galax y catalog ues created from the Millennium 
Run Simulation. |Paz et al. ( 2008. ) correlated galaxy and halo 
shape with large scale structure and also claimed evidence 
for a significant correlatio n in both th eir numerical models 
and in the SDSS (DR-6). lLand et all (2008) examined the 
evidence for a violation of large-scale statistical isotropy in 
the distribution of projected spin vectors of spiral galaxies, 
using data on line of sight spin direction classified by mem- 
bers of the public partaking in the online project "Galaxy 
Zoo". They found no evidence for overall preferred hand- 
edness of the Universe, even thou gh the related stu dy of 
spin correlations in G alaxy Zoo bv lSlosar et al.l (|2009l ) and 
Ijimenez et al.l (|2010l ) did find a hint of a small-scale align- 
ment^ 

iLee fc Erdogdu 1 1|2007^ determined tidal fields from the 
2MASS redshift survey (2MRS) and looked at the ori- 
entation of galaxies in the TuUy Galaxy Catalogue rela- 
tive to this tidal field. They found clear evidence of cor- 
relation between the galaxy spins and the intermediate 
principal axes of the tidal shear (see appendix [Aj. This 
ties in with orientation studies based on galaxy catalogues 
which indicate that the spin vector of galaxies tends to be 
aligned with the plane of the wall in which the halo is em- 
bedded jK ashikaw a fc Okamuralll992l : iNavarro et"al] |2004| : 
iTruiillo et al. 2006 ). 

Paradoxically, there are no systematic studies of the 
orientation of galaxies in filaments despite their relatively 
easier identification. We plan to use the edge-on galajcies in 
a filament catalogue in a marked correlation analysis. 



1.2.1 Spin alignment relative to structure 

Methods based on correlations between pairs of galaxies do 
not identify systemic large scale alignments or correlations 
of alignments with specific types of structures such as clus- 
ters or filaments: they merely tell us that the spin orienta- 
tions are not random. The main problem here is to accu- 
rately iden t ify st ructures and the objects that they contain. 
iHan et al.l (|l99d ) examined a cosmic filament of galaxies, 
but could find no evidence for any sort of alignment. 

Recent advances in cosmic web structure analysis have 
provided methods whereby galaxy redshift catalogues can 
objectively be segmented into voids, fil ament and clus t ers. A 
first step in this direction was taken bv lTruiillo et al.l (|2006l ) 
who, using the SDSS and 2dFRS catalogues, found evidence 
that spiral galaxies located on the shells of the largest cosmic 
voids have rotation axes that lie preferentially on the void 
surface. Their result tied in with numerica l expectations of 
such an alignment bv lBrunino et al.1 (|2007l ). 



1.2.3 Previous history 

The earliest commentaries on this subject (|Brownl (|l938l ) 
and Reynolds ( 1922) citing even earlier work by Brown) pre- 
date the establishment of an extragalactic distance scale by 
Hubble. Brown reported the detection of a significant sys- 
temic orientation of spiral galaxies having small inclinations 
to the line of sight. Much later, after the publication of the 
Palomar Sky S urvey , Brown publi shed two further papers 
on the subject (|Brownlll964l . 119681 ) based on visual exami- 
nation of thousands of galaxies covering v ast areas of sky , 
reporting evidence for significant alignment . iReavesl l|l958l ). 
using specially prepared plate material, was unable to con- 
firm Brown's claims. 

Brown's work was later taken up bv lHawlev fc Peebles! 

(jl975l ) who removed possible systemic biases by using a dou- 
ble blind experiment. They searched for evidence of over- 
all alignment, alignment between neighbouring pairs and 
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alignments in the Coma cluster of galaxies. None of these 
searches yielded statistically significant evidence for galaxy 
alignment. 

But that was far from the end of th e story. The Lo- 
cal Su p ercluster (LSC) wa s analy sed by Jaaniste fc Saar 
' 19781) ■ iMacGillivrav et al.l (|l982l i and iFlin fc Godlowski 



19861 ) with contradictory results, and with no significant 



evidence emerging. However, the latest attempts at resolv- 
ing the question o f align men ts in the Local Supe rcluster are 
those of lHu et al.l l|2006l) and lArval et al.1 (|2008al lbl) in which 
it is claimed that there is evidence of significant alignment 
of some galaxy sub-populations, but not others. The situa- 
tion from these surveys is manifestly unclear. The best that 
can be said is that, even if there is any systemic alignment, 
it is far from easy to find evidence for or to quantify the 
phenomenon using approaches based on analysing samples 
having large numbers of galaxies. 



1.3 A clean filament catalogue 

In this paper we try an alternative approach in which we 
build a small but statistically useful catalogue of filaments 
containing edge-on galaxies. The use of edge-on galaxies sim- 
plifies the problem of determining the space orientation of 
the galaxies. The filaments are places where, a priori, we 
might expect to find an alignment effect. Our catalogue is 
not in any way biased with respect to the orientations of the 
galaxies defining those filaments. 

The plan of the paper is to review the theoretical no- 
tions underlying spin orientation correlations with large 
scale structure, and comment on what N-Body experiments 
tell us in this regard. Then we go on to discuss how, using 
the DR5 release of the SDSS, we build a small catalogue 
of cosmic web filaments that are "good places to seek out 
orientation anomalies" . 



1.4 Outline 

We start by discussing the context of our study in section [S] 
the expected alignment of the angular momentum with the 
surrounding weblike structures. Subsequently, we present 
our SDSS DR5 galaxy sample in section [3l along with the 
applied corrections for artefacts such as redshift distortions. 
Section |3] elaborates on the construction of our FILCAT-0 
filament catalogue, and the selection of the class A filaments 
which will be subjected to our feature analysis. On the basis 
of the geometry of the galaxy-filament alignment configura- 
tion, in section [5] we will derive a statistical distribution 
function for the random galaxy-filament orientation distri- 
bution which forms the reference model for our study. The 
galaxy orientations of edge-on galaxies in our selection of 
SDSS filaments is analysed in section |6l followed by the sta- 
tistical feature measure analysis in section [T] Armed with a 
firmly established significant alignment of a small subset of 
14 galaxies, in section |8] we address the question of the na- 
ture of these aligned objects. Finally, section [5] summarizes 
and discusses our results. 



2 TIDAL MANIFESTATIONS: 



ALIGNMENTS AND THE COSMIC WEB 

Within the context of gravitational instability, the gravita- 
tional tidal forces establish an intricate relationship between 
some of the most prominent manifestations of the structure 
formation process (for a proper analytical definition of the 
tidal shear field Tij we refer to appendix . 



2.0.1 the Cosmic Web 

Perhaps the most prominent manifestation of the tidal shear 
forces is that of the distinct webUke geometry of the cos- 
mic matter distribution. The Megaparsec scale tidal shear 
forces are the main agent for the contraction of matter into 
the sheets and filaments which trace out the cosmic web. 
The main source of this tidal force field are the compact 
dense and massive cluster peaks. This is the reason be- 
hind the strong link between filaments and cluster peaks: 
filamen ts should b e seen as tidal bridges between cluster 
peak s l|Bond et aL _ 19961: van de Weveaert fc Bertschingeij 



1 19961 : Ivan de Weveaert fc Bondll2009l ) 



2.0.2 Angular Momentum and Tidal Torques 

In addition to the filamentary cosmic web itself, we recog- 
nize the manifestations of the large-scale tidal fields over 
a range of scales. Perhaps the most important influence is 
their role in generating the angular momentum and rotation 
of dark haloes and galaxies. It is now generally accepted that 
galaxies derive their angular momentum as a consequence of 
the action of tidal torques induced by the surrounding mat- 
ter distribution, either from neighbouring proto-galaxies or 
from the large scale structure in which they are embedded. 

As a result of the tidal force field, a collapsing halo 
will get torqued into a rotating object. The magnitude and 
direction of the resulting angular momentum vector Li of 
a protogalaxy is related to the inertia tensor, lij, of the 
torqued object and the driving tidal forces described by the 
tidal tensor Ti^ 



(1) 



with eijfc the Levi-Civita symbol, and where s ummation is 
implied over the repeated indices l|Whitell 19841 ). 

The idea that gala xy spin orig inated through tidal 
torques originated with iHovld (| 19491 ) who suggested that 
the source of tides was the cluster in which the galaxy was 
embedded. This was put into th e con text of modern struc - 
ture formation by Peeblei (|l969l ) and iDoroshkevicbl (|l970l ). 
where iPeeblesI l| 19691 ) focussed on the tidal interactions be- 
tween neighbouring proto-galaxies as being the source of the 
rotation. 

However, confusion concerning the efficiency of the 
mec hanism remained wid e spread until the numerical study 
by lEfstathiou fc Jones ] (|l979l ) and many generations 
of simulations thereafter (e.g. iBarnes fc Efstathio ul Il987l : 
iDubinski et"aLlll993l : iPorcianil [20021 : iKnebe fc Poweji2008l ). 
Although the early simulations did not attempt to iden- 
tify the source of the tidal field, numerical simulations 
along with related analytical studies demonstrated the vi- 
abilit y of the tidal torque mechanism ([Whitd 1984'; ^Rvdenl 
19871: iHeayens fc Peacoc^ 1 19881 : ICatelan fc Theu ns 199i; 



Lee fc Penll200ll : |Porcianil(2002l ) 
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2.0.3 Aligning Spin and Web 

The tidally induced rotation of galaxies implies a link be- 
tween the surrounding external matter and the galaxy for- 
mation process itself. With the cosmic web as a direct mani- 
festation of the large scale tidal field, we may therefore won- 
der whether we can detect a connection with the angular mo- 
m entum of galaxies or ga laxy halos. The theor etical studies 
of lSugerman et all l|2000l ) and lLee fc Peiil (|2000l ) were impor- 
tant in pointing out that this connection should be visible 
in the orientation of galaxy spins with the surrounding large 
scale structure. 

The notion that galaxy spins might be correlated with 
large scale structure has been extensively discussed on the 
basis of data catalogues, even those that existed in the early 
days of large scale studies (see sect. I1.2| l. Theoretically, the 
interpretation of this alignment is facili tated by invoking the 
parameterized formalism forwarded bv lLee fc Peril l|2000h to 
describe the correlation between the angular momentum L 
feqn. [Blj) and the tidal tensor Tij. To this end, they intro- 
duced the parameter c to quantify the autocorrelation tensor 
of the angular momentum in a given tidal field, 

(L,L,iT> cx i^.j + c{^5,, - T^kT^j) (2) 

where c — corresponds to the situation where tidal and 
inertia tensor would be mutually independent and the re- 
sulting angular momentum vectors would be randomly dis- 
tributed. In appendix[B]we summarize the main ingredients 
of this description. 

The alignment of galaxy haloes, the spins of 
galaxies and dark halos, of clusters and of voids 
with the surrounding large scale structures h ave been 
the s ubject of nume r ous st udies (see e.g. Binggelil 
19821: Imiee fc Katger tI Il987l: i Plionis fc Basilako^ 



Plionis et al. 12003: B asilakos et al . 2006; Truiillo ct al. 200 
Altav et alTT^Oel: lAragon-Calvo et al. 2007: IHahn et al 



'200d; 
200i 



Lee fc EvrardI l2007l : iLee et al., .2008 : .Platen et al. 
The expectation of a strong correlation between 
the orientation of structures and the large scale force field 
configuration is based on the existence of s uch correla- 
tions in the primordial G aussian density field (|Bondl [l987l : 
iDesiacgues fc Smithl200^ . Conversely, one should recognize 
the important role assigned to the alignment of rich clusters 
in de termining th e strength and morphology of the cosmic 
web jBond et al.|[r996i ). 



2.0.4 Aligning Spin and Web: 

evidence from N-body simulations 

N-Body simulations provide information on the alignment of 
the orientations of galaxy halos and galaxy spins. With this 
it is possible to discuss the correlations of spins among galax- 
ies and the correlations of spin with large scale structures 
such as filaments and walls. The main difficulty in the latter 
case is to unambiguously identify those large scale struc- 
tures. The recent advances in cosmic web structure analysis 
have enabled a systematic segmentation of the cosmic mat- 
ter distribution into voids, sheets, filament and clusters. 
As a resul t, r ecent A^-body simulations 



sheetlike nature of the environment has a distinct infiuence 
on the shape and spin orientation of dark matter haloes. 
In the case of of haloes located in large scale walls, they 
seem to agree that both the spin vector and the major 
axis of inertia lie in the plane of the wall. In the case of 
the alignment of halos w ith t heir embedding fila ments, 
lAragon-Calvo et al.l |20o3) and iHahn et all (|2007al ) found 
evidence for a mass and redshift dependence , which has 
been confirm e d by the studies of |Paz et al] (|2008l ) and 
IZhang et al. I (I2OO9I ). This mass segregation involves a 
parallel alignment of spin vector with the filament if the 
mass of the halo is less than the characteristic halo mass 
m < 10^^ M0 h~^, turning into a perpendicular orientatio n 
for more massive haloes (also see iBailin fc Steinmet j 2004 ) . 
A tan talizing related finding is that bv iFaltenbacher et al] 
(|2009t ). who claims correlations between spin and large 
scale structure extending over scales of over 100/i"^Mpc. 

Results from A^-body simulations can not be directly 
compared with observations. The reason for this is that the 
only available information of the spin of galaxies comes from 
the luminous baryonic matter. The coupling between the 
spin of the baryonic and dark matter components is still 
not well understood. One of the key issues raised by studies 
of simulations is that the dark matter halos and the bary- 
onic disk component may not line up with one another. 
Recent high-resolution hydrodynamical simulations of the 
forma ti on of galaxi e s wit hin a cosmic filament by iHahnl 
(|2009l) : iHahn et al.1 l|2010h co nfirms, or even stren g thens , 
the earlier trend obse rved by Ivan den Bosch et al.1 (I2OO2I ) 
and IChen et al.1 (|2003D . who found that a median angle be- 
tween gas and d ark matter is in the order of ~ 30° (also see 
iBett et al.1 l2009il This misalignment between dark matter 
and gas is sufficient to erase any primordial alignment signal 
between the galaxies and their host filament or wall. This 
effect however, could be rendered negligible if the gaseous 
component of galaxies retains its primordial orientation bet- 
ter than its dark matte r counterpart as sug gested by the 
SPH simulati on study of Navarro et al.l ||2004|). The fi lament 
simulation bv lHahrJ ( 2003 ) (also see lHahn et al.ll2O10h seems 
to be considerably less optimistic in this respect. 



3 THE SDSS GALAXY SAMPLE 

The results presented in this letter are based on a galaxy 
sample selected from the largest galaxy survey to date, the 
Sloan Digital Sky Survey data release 5 (SDSS DR5). The 
SDSS is a wide-field photometric and spectroscopic sur- 
vey carried out with a de dicated 2.5 meter telescope at 
Apache Point, New Mexico (lYork et al.| . 20001. The telescope 
scans continuously the sky on five photometric band-passes 
namely u, q, r, i and z, down to a limiting r-ban d magnitude 
of 22.5 (|Fukugita et al.lll996l : [Smith et al.ll2002l ). 

The analysis presented in this letter includes approxi- 
mately 100, 000 SDSS DR5 galaxies located in the range 

< ?7 < 40 
0.01 < 2 < 0.11 



(lAragon -Calvo et al.l |2007|; [H ahn et_ al.l l2007allbl: iPaz et al.1 
[2OO8; H ahn et al.ll2009l : lHah n 2009; Z hang et all2009l ') have 
been able to find, amongst others, that the filamentary or 



within the survey coordinate system (A, r?) (jStoughton et al.l 
I2OO2I ). This geometry encloses a volume of approximately 
7.3 X lO'^ Mpc^h-^ 

In order to take full advantage of the galaxies in the 
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Figure 1. The SDSS DR5 sample volume seen in two projections. Top: redshift slice through the volume, out to 2: = 0.1. The image shows 
the SDSS galaxies in the sample volume that have been identified by MMF as filament galaxies. The galaxies are shown at the location 
assigned by the MMF filament compression algorithm (see Aragon-Calvo et al. 2007). The small open circles indicate the location of all 
492 edge-on galaxies in our sample. The 14 large ellipses, and corresponding number, locates the 14 edge-on galaxies that our analysis 
identifies as significantly aligned with their parent filament. The elliptical shape corresponds to the orientation of a circle with a similar 
orientation as the galaxy. Bottom: sky region of the SDSS DR5 survey, including the location of the 14 positively aligned galaxies (see 
fig. 19 for galaxy images). The edge-on galaxies are plotted as lines indicating the galaxies' sky orientation. Each of these galaxies is 
shown against the backdrop of the projection of their embedding filament. 



sample we used a magnitude limited catalogue. This means 
that the radial distribution of galaxies must be weighted in 
order to produce an isotropic distribution. We used a simple 
formula to model the change in the me an number of galajcies 
as fun ction of their redshift given by lEfstathiou fc Moodvl 
(|200lf ): 

dN = Az^exp{~{z/zrf). (3) 



Where A is a normalization factor that depends on the den- 
sity of galaxies, Zr is the characteristic redshift of the distri- 
bution and /3 encodes the slope of the curve. 

Before applying filament finding algorithms, like MMF, 
it is necessary to fill in any holes in the survey (see sec- 
tion I3.1.2P and handle the effects of the boundaries of the 
surveyed region (see section [3.1.3p . 
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3.1 Data systematics, removal of artefacts and 
the final filament samples 

The catalogued galaxy distribution is affected by several 
systematics and artefacts. Among the most important are 
telescope artefacts, unsampled areas of the sky, peculiar ve- 
locities from both linear and non-linear processes. 



is extreme). Using the axial ratio criterion we can measure 
the inclination of face-on and edge-on galaxies with an un- 
certainty of ±12°. Imposing a stricter selection criterion on 
the axial ratio would improve this, but it would also decrease 
the size of the useful sample of galaxies. 

Within the boundaries of our SDSS DR5 subsample we 
found 492 edge-on galaxies. 



3.1.1 Fingers of God 



Fingers of God have a characteristic elongated shape that 
will introduce "false" filaments in the distribution of galax- 
ies along the line of sight. In order to correct this effect 
we compressed the finge rs of God in a similar way as done 
bv lTegmark et ahl (|2004 ) by identifying elongated structures 
and them compressing them along the line of sight. One 
drawback of the finger of God compression algorithm is that 
it systematically removes filaments oriented along the line 
of sight. This is a consequence of those filaments being con- 
fused with fingers of God. In our analysis we take account 
of this. 



3.1.2 Unsampled regions 

We correct the holes in the angular mask by means of a new 
method that exploits the volume fi lling properties of the 
Delaunay tessellation (|Arag6n-Calvo 2007). The Delaunay 
Tessellation Field Interpolator is a geometry-based interpo- 
lation scheme that uses the spatial information at the edges 
of the hole and follows the tessellation to derive the density 
field in the interior of the hole. 



3.1.3 Boundaries 

In order to produce an uniform density field around the sur- 
vey volume we place particles, randomly, with the same ra- 
dial distribution as the galaxies in the sample, outside the 
sample volume. 



3.1.4 Edge-on selection 

Optical surveys such as the SDSS provide enough informa- 
tion to derive the position angle, inclination and, in some 
cases, the sense of rotation of the galaxy. However we can- 
not determine which is the approaching side or the receding 
side. This results in a four- fold degeneracy in the derived 
spin vector. The study of spin alignments only requires the 
direction of the spin vector and not its sense of rotation. This 
reduces the degeneracy to two-fold which is still insufficient. 

There are two special cases where the spin vector can 
be unambiguously determined only from the position angle 
and inclination: i) Edge-on galaxies having their spin vector 
in the plane of the sky. ii) Face-on galaxies having their spin 
vector pointing along the line-of-sight. 

In this paper we restrict our analysis to edge-on galax- 
ies. Edge-on galaxies are identified in terns of the ratio of 
their projected axes by r^/ra < 0.2. Strictly speaking, this 
applies only to spiral galaxies which are assumed to be fiat 
rotating disks with their spin vectors pointing perpendicu- 
lar to the disk. (We could use colour information to separate 
early and late type galaxies if that were necessary, but we do 
not do that since the fiattening of the objects in our sample 



4 FILAMENT CATALOGUE: 

DETECTION AND SELECTION 

Even if there were systemic, non-random, orientations set up 
by, say, global tidal fields at the time of galaxy formation, 
subsequent dynamical evolution would conspire to largely 
eradicate the phenomenon from the present data. Here we 
think of interactions between galaxies among themselves or 
with neighbouring large scale structures such as clusters. 
Our approach, then, is to look for fossils that have not un- 
dergone such evolution and so retain a vestige of their orig- 
inal orientation. 

To accomplish this we use techniques that are perhaps 
more familiar in data mining: we cull and clean the dataset 
to obtain a subset of the very best data that is available, and 
then analyse the result to look for potentially interesting 
anomalies in the distribution of orientations. We address 
the general significance of context of this strategy below, in 
sect. WT\ 

We use a scale free structure finding algorithm (MMF, 
see below) to identify cosmic web filaments, restricting our- 
selves to the edge-on galaxies within those filaments (see 
above) . The galaxy spin axis orientation is then in the plane 
of the sky, and since MMF yields the space orientation of 
the host filament, we can calculate for the entire sample the 
distribution of spin axis - filament angles. The goal is then 
to answer the question as to whether these orientations are 
random, and if not, say whether there are preferred orienta- 
tions. 



4.1 Sample analysis: KDD data mining 

The current approach differs substantially from previous at- 
tempts to find evidence for alignments between galaxy spins 
and large scale structures in that we search, not for a trend 
or correlation in the data, but for specific candidate exam- 
ples of such alignments. By focussing on edge-on galaxies 
in filaments, and then grading the resultant filaments we 
arrive at a relatively small database of high quality fila- 
ments hosting edge-on galaxies that can be analysed from 
the point of view of relative galaxy-filament orientation. To 
do this we follow the precepts of the so-called KDD process 
for analysing large, possibly heterogeneous or partially sam- 
pled, data sets. It involves several steps including culling and 
cleaning of the search data set, in this case the DR5 release 
of SDSS. 

KDD stands for "Knowledge Discovery in Databases" , 
a term intr oduced by Piatetskv-Shapirol l|l99ll ) . KDD is de- 
scribed bv lFavvad et al.l l|l996l )" as "The non-trivial process 
of identifying valid, novel, potentially useful, and ultimately 
understandable patterns in data". The more familiar con- 
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Figure 2. A selection of 6 filaments showing the compressed filament and the location of the host galaxy. Each panel shows the 
identification of the filament in the FILCAT-0 filament list and the assigned classification. Black dots: galaxy location after filament 
compression (Aragon-Calvo 2007). Grey dots: original SDSS location. 



cept of Data Mining is very much a part of KDD, rather 
than the other way arounqj. 

An integral part of the KDD approach is to have a mea- 
sure of "interestingness" that ranks what we are looking for. 
Note that this is not necessarily the same as the more fa- 
miliar concept of significance. In this respect we study the 
histogram of the distribution of spin orientation of a galaxy 
relative to the parent structure, identifying features in the 
histogram that make it notable, or interesting, with respect 
to the problem at hand. We introduce a "Feature measure" 
that reflects this level of interestingness. 

Formally, the process of extracting a suitable feature 
measure for the present study is based on wavelet analysis 
of the histogram of the distributions of orientations from the 
data sample itself and from a simple simulation. 



Web. Several methods have been used in an attempt to 
identify and extract the morphological components of the 
Megaparsec-scale m atter dis tribution ijBarrow et al.l Il985l: 
Babul fc Starkmanl[ l992: Lu o fc Vishniadll995l:IStoica et all 



20051 : IColberg et al.ll2005i : iPimbletll2005riSousbie et al.ll2008l . 



among others) with varying degrees of success. 

Key to our approach is to have a first class structure 
finder so that we can unambiguously identify the filaments 
and the member galaxies, and so find any edge-on galax- 
ies they might be hosting. It should be stressed that in any 
study of galaxy orientation relative to parent structure, it is 
important to identify that parent structure as accurately 
as possible. The results presented here are based on the 
Multiscale Morphology Filter MMF (,Arag6n-Calvo et all 
Based on a densit y field reconstructed from SDSS 
using the DTFE process Jschaap fc van" de Wevgaerj|2000l : 
Ivan de Weveaert fc Schaap 1 12009|) . MMF largely achieves 
the goal of outlining the filamentary spine of the Cosmic 
Web. However, as we shall see, it is worth culling the fila- 
ment sample to provide a filament sample that is most suited 
to this task. We want to avoid dilution of what, in the light 
of past studies, is evidently at best a weak efi^ect. 

The morphological characterisation achieved by MMF 
enables us to isolate specific host morphological environ- 
ments for galaxies and test predictions from the tidal torque 
|http : //www . cs ■ csl ■ cuny . edu/ - Imberman/DataMlnlng/KDDbeginnlngs . ^ibSj Dry in a systematic way. 



4.2 



MMF Filament identification and 
classification 



Revealing correlations between galaxy spin direction and 
the large scale structure requires the ability to unambigu- 
ously identify the morphological features of the Cosmic 



^ There is a good lecture by Susan Imberman at 
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Figure 3. Specific example of an edge-on sample galaxy, nr. 347 of our sample. Bottom right; SDSS image of edge-on sample galaxy nr. 
347. Top right: location of galaxy nr. 347 on the SDSS DR5 sky region and redshift region, along with some key data. The position and 
orientation of the galaxy with respect to its embedding filament is depicted by means of three mutually perpendicular squares of 20h~^ 
Mpc size, centered on the galaxy. Black dots; galaxy location after filament compression (Aragon-Calvo 2007). Grey dots; original SDSS 
location. 



4,2.1 The Multiscale Morphology Filter 

The MMF objectively segments a point set representing 
the cosmic web into its three basic morphological compo- 
nents: clusters, planes and filaments. It does this first by us- 
ing the DTFE methodology to create a continuous density 
field from the point set, and then hierarchically analysing 
the properties of the local matter distribution. The DTFE 
method produces an optimal reconstruction of the contin- 
uous density field, retaining the characteristic hierarchical 
and anisotropic nature of the Cosmic Web. This allows 
MMF, multi-scale by construction, to produce a catalogue 
of filaments with relatively few scale biases. 

The Multiscale Morphology Filter is based on the 
second-order local variations of the density field as encoded 
in the Hessian matrix (d^p/dxidxj). For a given set of 
smoothing scales we compute the eigenvalues of the Hessian 
matrix at each position on the density field. We use a set of 
morphology filters based on relations between the eigenval- 
ues in order to get a measure of local spherical symmetry, 
filamentariness or planarity. The morphological segmenta- 
tion is performed in order of increasing degrees of freedom 
in the eigenvalues for each morphology (i.e. blobs — >■ fila- 
ments). 

The response from the morphology filters computed at 



all scales is integrated into a single multi-scale response 
which encodes the morphological information present in the 
density field. In a filament, the eigenvectors of the Hessian 
matrix corresponding to the smallest eigenvalue ep) indicate 
the direction of filament. Eigenvectors are computed from a 
smoothed version of the density field in order to avoid small- 
scale variations in the direction assigned to filaments. 

A little more detail is presented in Appendix: ICl for the 
sake of completeness. 

4.2.2 The raw FIL CAT- filament catalogue 

On the basis of our MMF technology we construct a raw 
filament catalogue. Before this can be used it is necessary 
to correct for the effects of redshift distortions and to go 
through a filament thinning (or compression) process to en- 
hance the visual appearance for quality assessment (see sec- 
tions 13.1. 111 . The final stage is to grade filaments and se- 
lect those deemed most useful for angular momentum fossil 
searches (see section 1^3)) . 

The MMF method yields a sample of 426 filaments con- 
taining edge-on galaxies that we shall use in the analysis 
of the spin-filament relationship. We refer to this catalogu e 
as the FILCAT-0 filament catalogue (| Aragon-Calvo II2007I '). 
Each filament can be viewed in three orthogonal projections, 
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and for each filament there is a raw view (redshift space) and 
a processed view in which a correction for redshift distortion 
has been apphed, followed by a compression of the member 
galaxies towards the spine of the filament. 

Provided the selection takes no account of the orienta- 
tion of the edge-on galaxy this procedure should be relatively 
benign. Accordingly, the classification was made in such a 
way that it is independent of the orientation of the sample 
galaxies, and also does not involve prior knowledge of the 
orientation of the filaments with respect to the line of sight. 



4.3 Filament classification 

Since MMF simply classifies structure on the basis of shape, 
the FILCAT-0 catalogue inevitably contains a variety of 
structures, all having a general linear structure. In the foren- 
sic spirit of this investigation we have to pare this list down 
to a subset of filaments where we can reasonably expect 
to find evidence for any putative relationship between the 
edge-on galaxy and its host filament. 

Filaments were scored 1, 2, or 3 depending on the visual 
appearance of the filament on three orthogonal projections 
of the redshift distortion corrected and compressed filament. 
Scoring was done independently by BJ and RvdW, and the 
scores totalled. Filaments scoring '1' by both, i.e. having a 
total score of 2, were assigned to "Class A", those with a 
total score of '3' to "Class B", those with score '4' to "Class 
C" and the rest to "Class D" . 

The 1-2-3 classification of filaments depended on a sub- 
jective quality assessment of the filament and on the location 
of the edge-on galaxy relative to the filament. The following 
criteria were used for filament quality: 

(i) must be defined by a reasonable number of galaxies 

(ii) must not be too sparse 

(iii) must not have a branched structure 

(iv) must be clearly defined in all three projections 

This refiects the use to which the filament is to be put: we 
are looking for places where we might reasonably expect 
the local tidal dynamics to be simple. In keeping with that 
approach we further require that the key edge-on galaxy: 

(i) must be located on or near the filament 

(ii) must not lie in a sparse area of the filament 

(iii) must not lie at the end of the filament 

The procedure in applying these criteria is as follows. For a 
"1" score the filament must be clearly defined in all three 
planes, and it must have a good population of points. The 
logic is that if it is not seen in all three planes, it might be 
a wall. Thus if this criterion fails in one plane it is ranked 
as a "3" . If the filament is well defined but sparse it is given 
a "2". Finally "1" and "2" filaments are downgrade to a 
"3" if the galaxy is not on a filament (subjective!) or if it 
sitting at the end or junction of filaments (a Y-shape). This 
is not unlike the old-fashioned 1950's and 1960's approach 
to classifying objects like clusters of galaxies. 

Although the orientation relative to the line of sight 
was not a classification factor, it turned out that all fila- 
ments making an angle closer than about 40° with the line 
of sight were eliminated. This is shown in figure [5] where 
the distribution of angles the filaments make with the line 




Figure 4. Vectors and angles defining the orientation of the fila- 
ment to the Une of sight and the orientation of the galaxy minor 
axis on the plane of the sky. The line of sight looks along the 
_ftr-axis and the galaxy minor axis defines the /-axis in the p[lane 
of the sky. The vector J completes the axis triad relative to which 
the filament orientation is measured. 

of sight is shown for the entire sample and for the Class A 
sub-sample. 

The final list of Class A filaments contained 69 objects 
only from the original sample of 492 edge-on galaxies. There 
are 106 Class A and Class B filaments, containing in total 
121 edge-on galaxies. 

5 GEOMETRY OF FILAMENT AND GALAXY 
ORIENTATIONS 

5.1 Angles 

Consider one edge-on galaxy in our sample. The spin axis 
of the galaxy (taken to be the minor axis of its shape) and 
the line of sight define a right handed coordinate system as 
shown in figure The filament will make an angle rj with 
the line of sight. The projected angle between the galaxy 
spin axis and the projection of the filament on the plane of 
the sky will be denoted by $. The angle of interest for this 
study, 9, is the angle between the filament and the galaxy 
spin axis. 

There is a four-fold degeneracy in the definition of these 
angles. Since the galaxy is only viewed edge on the spin 
axis S lies in the plane of the sky. There is no additional 
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information about the sense of rotation. The spin axis is 
defined without a sense of direction and so we should be able 
to make the transformation S — ^ — S without changing the 
angle 9. The line of the filament, F is also defined without 
a sense of direction and so making the transformation F 
—F should not change the value assigned to 6. 

We can handle this degeneracy by restricting the ranges 
of the angles to ^ ^ 7r/2 and ^ $ ^ tt, and by forcing 
the range of S to be ^ 7r/2. 

The angle 6 between the filament and the galaxy spin 
axis is defined in terms of the observables r] and $ by 



cos 6 — I sin 1] cos $1 



(4) 



The absolute value is needed because of the degeneracies 
discussed above and to bring 6 into the range ^ 6 ^ n /2. 

It is the statistical distribution of the angle 6 that con- 
cerns us. We are comparing the distribution of this angle, or 
specifically the distribution of cos(S), with what would be 
expected if the spin axis and the spine of the filament di- 
rection were uncorrelated. So, for each galaxy and its local 
filament direction, we compute cos(S). 



5.2 Uncorrelated spins and filaments 

The angles and <I> are clearly independent. In the absence 
of selection effects, cosr; is uniformly distributed, as is the 
angle $. Note, however, that if there were a correlation be- 
tween the directions of the spin axis and the filament, $ 
would no longer be uniformly distributed. If the spin direc- 
tions of galaxies are independent of the host filament, and 
the host filaments are themselves randomly oriented with 
respect to the line of sight, then the distribution of cos^ 
will also be uniform. 

However, owing to a variety of issues, the filament sam- 
ple is not randomly oriented and so cos?7 is not uniformly 
distributed (see figure[5]). Consequently, cos 9 will not be uni- 
formly distributed either. We see in figure [5] (bottom) that 
subsamples A and B of filaments have an almost uniform 
distribution of cos 77, except for a cutoff' at higher values of 
COS?;. For our sample, this corresponds to an exclusion of all 
filaments within an angle of rjc ~ 35° — 40° to the line of 
sight (see sect. 16.11) . 

If this is modeled with a sharp cutoff at angles -q < r]c, 
the distribution of cos 6 can easily be derived in closed form. 
The expected fraction of objects having a particular value 
of cos^, relative to a uniform distribution is 



/ — — arcsin 

TT 



cos 77c 



sin 6 



= f 



6 > 1x12- rj^ 
otherwise. 



(5) 



This is derived and compared with simulations in Appendix 

El 

Other factors come into assessing the distribution that 
would be seen even in the absence of a spin-filament corre- 
lation. In particular, the sample of filaments found by MMF 
with the corrections for redshift distortion and the compres- 
sion procedure produces a non-uniform distribution of cos(77) 
that is biased against including filaments lying along the line 
of sight. 

Subsequent visual selection of a sub-sample of filaments 
imposes additional constraints on the expected distribution. 
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Figure 5. Distribution of the angle, r], the filaments in our sample 
make with the line of sight. Bottom: the cost; distribution of 
the 70 class A filament galaxies (solid red bars) and i21 class B 
filament galaxies (solid blue bars), along with the distribution of 
the edge-on filament galaxies embedded in these filaments (dashed 
bars). Top: the distribution of the cos 77 values of the embedding 
filaments of the entire sample sample of 492 edge-on filament 
galaxies is shown in grey along with three sub-samples: the 70 
class A filament galaxies (dark blue) and the set of 121 Class B 
filament galaxies (red) and all edge-on galaxies in filaments not 
rated as '3' by either classifier (yellow). 



As it turns out, these can be modelled and so we can com- 
pare the data with a well-reasoned model. 



6 DATA ANALYSIS 

We can derive for our sample of filaments, and for a variety 
of sub-samples, the angle that the filament makes with the 
line of sight. This is shown in figure (5] 
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Figure 6. Distribution of the angle 6 the spin axis of the galaxy makes with the direction of its parent filament. Two distributions are 
shown: one the distribution expected from the analytic result of equation (??) and the other for the galaxies in the sample of Class A 
filaments. The hump in the filament data at cos 9 < 0.2 is outstanding. 

(top panel). The drop-off at cost; > 0.8 is hardly surpris- 
ing in view of the overall distribution: it presumably reflects 
the prejudice against finding filaments lying along the line 
of sight. 

The distribution of cos 77 for the sub-samples of Class 
A and Class A-l-B filaments is uniform for cos 77 < 0.8, as 
would be expected for a randomly oriented vector relative to 
the pole of the distribution (see lower panel, fig. O solid red 
and blue bars). So, while the distribution of the line of sight 
angles for the whole sample is far from uniform, the Class 
A and selection criteria have led to a uniform distribution 
of angles that has been censored at angles less than around 
35° — 40° to the line of sight. Not only does this finding 
assure us that we understand the angular distribution of 
filaments (see sec.|D}, but also provides a strong argument 
for the solidity of the MMF filament finding process. 

Because of this distribution of rj, the distribution of 
cos6, the angle the spin axis makes with the filament, is 
not uniform either. This readily follows from the geometric 
model described by equation [S] 

6.1 Filament angles with line of sight 

6.2 Orientations in our filament sample 

First we should remark that the orientations of the spin axes 
relative to the sky are consistent with a uniform distribution: 
there is no bias towards specific orientations in the SDSS 
catalogue. 
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Figure 7. Relationship between the inclination of the filament 
to the line of sight, jj, and the relative orientation of the galaxy 
and the host filament 9 for the Class A sample. The horizontal 
lines delineate the histogram bins shown in figure [6] The excess 
at low cos d is seen to be associated with filaments having sin rj > 
0.75, i.e. tending to be transverse to the line of sight. This should 
occasion no surprise as most of the Class A filaments turn out to 
be so-orientated. 



In figure [5] we see the decline in the number of fila- 
ments as cos ?7 — >■ 1: the distribution for the entire sample is 
nowhere flat. It is clearly reflected in the angular distribu- 
tion of parent filaments of the 492 edge-on filament galaxies 
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The distribution of spin axis to filament angles, 9, • 
rived from our Class A sample is shown in figure [G] We 
a striking hump for objects with cosO < 0.2: this conta 
14 of the 82 galaxies while in an uncorrelated distribut 
we would have expected less than 1. 

It is on this basis that we draw two conclusions: 

(i) The spin axes of these edge-on galaxies are correlal 
with the direction of their local filament in the sense tl 
their spin axes are perpendicular to the filament spine. 

(ii) The filaments are of real physical significance (oth 
wise we would not see anything other than a uniform disi 
bution, as per the non-selected data sample). 

It is important to take account of the fact that there 
no filaments in our Class A sample making a small anj 
7] < 7r/4, with the line of sight (see figure [5]). Remov 
these from the sample further decreases the occupancy of 
the bins having spin-filament angles greater than 7r/4 (i.e. 
cosd < 0.7) (see equation ((Sjl). The distribution of cos 9 for 
our Class A filament sample is shown in figure |6l and it is 
compared with the expectations of our simple uncorrelated 
spins distribution in figures [TOl and [TT] 



6.3 The significance of the anomaly 

The cos6' < 0.2 feature in the distribution is quite strik- 
ing. Simply considering the geometry as depicted in fig- 
ure |4] it is difficult to imagine any model based on random 
spin-filament orientation in which the number of objects in- 
creased with decreasing cos 9: the available phase space for 
the angles rj and "1> does not allow that. 

There are only 14 objects in this hump in the distribu- 
tion, and so, having found it, the question arises as to why 
the other galaxies in the edge-on sample are not correlated 
in this way. Why are these 14 galaxies so different from the 
rest ? 

Before addressing this we must however reassure our- 
selves that we have done nothing to the data that might 
give rise to this effect. After all, we have created a very spe- 
cific sub-sample from some 426 filaments hosting edge-on 
galaxies. 

Several factors have come into the sample definition: 

(i) Corrections for fingers of god eliminate filaments along 
the line of sight. Since for such a filament the spin axis will 
tend to be perpendicular to the filament, this elimination 
actually removes objects having small cos 9 from the sample! 
So this does not help enhance our bump 

(ii) We have rejected filaments in which the galaxy in 
question does not lie close to the spine of the filament. In- 
cluding such galaxies would only serve to randomise the dis- 
tribution rather than enhance a feature at some particular 
angle. 

(iii) We have rejected filaments for which a spine cannot 
unambiguously be defined. This includes rejecting filaments 
that may be branching in or around the position of the edge- 
on galaxy. 



6.4 Other samples - a comparison 

We shall, in the next section, make a statistical assessment 
of this result by comparing it with simulated samples. It 
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Figure 8. Distribution of the orientation of the spin axis of an 
edge-on galaxy relative to the direction of the host filament, for 
the entire sample of filaments. 



is useful in the first instance to look at a variety of other 
samples. Firstly we shall randomise the spin vectors among 
the filaments in the Class A sample to show that this result 
disappears. Then we shall comment on whether the entire 
catalogue shows any evidence of this, finally we shall look 
at the edge-on galaxies in filaments that appear to lie in the 
plane of the sky. 



6.4-1 The entire sample 

When classifying the filaments we rejected around 70% of 
the original sample. It is interesting to compare the selected 
Class A sub-sample with the original sample from which it 
was extracted, remembering of course that the original sam- 
ple had a non-uniform distribution of filament orientations, 
and that the position of the edge-on galaxy was not con- 
strained. 

If we repeat the analysis for a selection of galaxies from 
poorly defined filaments we see a distribution that is con- 
sistent with the spins being uncorrelated. This is shown in 
figure [8] The distribution is almost uniform, but there is 
nonetheless a discernible feature at low cos{9) which we 
might convince ourselves is the same hump as in the se- 
lected sub-sample. However, unlike in the selected sample, 
this hump does not have a striking level of significance and 
by itself would not constitute particularly compelling evi- 
dence for alignment. 



6.4-2 Filaments in the plane of the sky 

Filaments in the plane of the sky (or close to it) offer another 
opportunity to test for alignments within this sample. In 
this case sin 77 « 1 in equation [3] and the angle between the 
filament and the spin axis is precisely what is seen. For such 
filaments the redshift distortion is presumably irrelevant, 
causing a broadening of the galaxy distribution. 

In order to maintain the quality of the sample we re- 
strict ourselves to sky-plane filaments rated '1' or '2' by both 
BJ and RvdW: i.e. we reject any filament rated '3' by either. 
This leaves a sample of 51 filaments. The distribution of cos 9 
for this sample is shown in figure [9] The sample is selected 
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Figure 9. Filaments in the plane of the sky; distribution of the 
angle 9 the spin axis of the edge-on galaxies make with the host 
filaments. The blue histogram is the complete sample of sky-plane 
SDSS filaments and the red histogram is a sub-sample excluding 
poorly rated filaments (see text). 
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Figure 10. Red: distribution of the angle, 9 in bins of equal 
cos 9 for the Class A sample of spin vectors. White: The expected 
distribution of spin axis directions when edge-on galaxies in our 
sample are assigned to randomised host filaments. 
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Figure 11. Red: distribution of the angle, 9 in bins of equal cos 9 
for the Class A sample of spin vectors. Grey: distribution 9 in 
bins of equal cos 9 expected according to a model for the filament 
selection process consistent with figure O 



primarily on the basis of the angle 77. It is perhaps surprising 
that there are only three filaments from the Class A sample 
having cos 9 < 0.2 in this sample of sky-plane filaments: the 
samples are effectively independent. The fact that we see in 
the sky-plane filaments a hump at low cos 9 is support for 
the conclusion suggested on the basis of the more rigorously 
selected Class A sample. 



7 STATISTICAL MODELLING 

The Class A sample consists of 67 filaments, of which we sin- 
gle out 14 for special attention. If we extend this to include 
the Class B filaments the sample grows to 121 filaments. 
Generating a simple model for random spin orientations is 
straightforward, but the question remains how to evaluate 
the credibility of the observed distribution. The "natural" 
test would be to compare the data with the simulation us- 
ing a test like the Kolmogorov-Smirnov (K-S) test. How- 
ever, such tests are generally sensitive to any differences in 
the distributions and have a low level of discrimination in 
indicating precisely what that difference is. 

What drew our attention to this result was the unex- 
pected peak of the spin orientation for values of cos6 < 0.2. 
Had we seen a peak in the range 0.2 < cos6 < 0.4 we might 
have been equally impressed. So the idea is to assess the 
observed distribution relative to a set of distributions that 
would have caught our attention. We do this by introducing 
"feature measures" that measure the features that single out 
the histogram as being of particular interest, and ask how 
often we might have seen one that was at least as interest- 
ing in this respect. The feature measure is perhaps more a 
measure of how justified we should feel in regarding the dis- 
tribution as perceptually exceptional, rather than a measure 
of statistical significance. 



7.1 Feature measures 

We calculate, for the binned histogram of the distribution 
of spin-filament angles, three indices based on the cos 9 his- 
togram (since it was the binned histogram that was the ob- 
ject that caught our appreciation). The feature measures are 



FMl 



FM2 
FM3 



6m(0.1) -f 6m(0.2) 

-2[6m(0.3) + 6m(0.4)] + 
-|-6m(0.5) -I- 6m(0.6) 
6m(0.8) + 6m(0.9) -f- 6m(1.0) 
|FA/1| +FM2 



(6) 
(7) 
(8) 



where bin{n) refers to the occupancy of the bin n — 0.1 < 
cos 6* < n with < n < 1.0. The use of these particular 
measures in describing what we see is rather self-evident, 
but they can be formally derived from the data using com- 
ponent wavelet analysis of the histograms of the simulated 
distributions and the data. 

FMl will be recognised as an estimator of the curva- 
ture of the interval < cos 6' < 0.6. FM2 is the occupancy 
of the last three bins. FAI3 is a composite index refiect- 
ing both these attributes. Note that we have constructed 
FM3 using the modulus of FMl: we want FM3 to refelct 
any histigram that might have aroused our attention, not 
only the histogram we have derived from the data. Using 
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|_FM1| rather than FMl in defining FM3 simply reduces 
the measured level of "interestingness" of our histogram. 
Since these feature measures are in fact rather broad aver- 
ages: they avoid sensitivity to details within the distribution 
and so they reflect broad features that would attract atten- 
tion when viewing data. 

The feature measure gives lower signiflcance levels than 
would a simple test like the K-S test since it allows a greater 
number of accidental possibilities to be considered as being 
similar to or greater than deviations from the reference sam- 
ple. 

We shall make two comparisons. The first with a sample 
using our SDSS Class A data in which the spin axis orienta- 
tions have been randomly reassigned to the filaments in the 
same sample, The second using a sample constructed from 
a model in which galaxies have randomly oriented spin axes 
relative to a sample of filaments having the same selection 
function as the observed SDSS sample. 



7.1.1 A shuffled Class A sample 

If we randomize the data by simply shuffling the spin axes 
in the sample among the filaments we get the distribution 
shown in figure 1101 The distribution of cos 6 is consistent 
with a uniform distribution. 

This would lead us to believe that there is indeed a 
strong effect in the data. For reference we also show the 
histogram that would arise if the spin axis orientation reltive 
to the host filament were random. We have chosen a flat 
distribution for cos rj with a cut-off at = 45° , as suggested 
by figure [5l 



7.1.2 A simulated sample 

We simulated a random distribution of spin axis - filament 
angles, 9, for a sample of filaments whose angle with the line 
of sight, 77, followed a simple distribution like that indicated 
in figure [5] for the distribution of our Class A filaments. We 
modelled this with a uniform distribution cut off at cos 77 > 
45°. 

From this distribution we randomly generated 1000 
samples, each containing 70 points with a value for 6, and 
calculated the feature measures FM1,FM2,FM3 as per 
equations ([B]), (O, Since FM2 can be positive or nega- 
tive, we used both FM2 and |_FM2| as measures. The first 
of these is indicator of a valley in the low cos 6 distribution, 
while the second of these looks for any signiflcant feature 
there, be it a valley or a hump. The feature measures for 
the random sample were then ranked: the position of our 
sample in these rankings was taken as the measure of signif- 
icance. 



7.1.3 First Remarks 

It will be noticed that the mean samples depicted in flgures 
1101 and rrn are not quite the same: the distribution shown in 
flgure [To] is somewhat flatter than that shown in figure 1111 
The former is based on a large number of random shuffles 
of the actual data while the latter is based simply on our 
model. The difference almost certainly lies in the simplicity 
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Figure 12. The absolute value of the curvature of the low cos 8 
bins as indicated by our index |FM1|. Our SDSS sample is rep- 
resented by the black spike, its height has no relevance except for 
clarity of display. Our sample would be ranked joint 8th (along 
with 4 others) relative to the total sample of 1000 cos 6 values: the 
two-sided probability of getting a feature attracting our attention, 
like that in our analysis, by chance is about 1% 
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Figure 13. The distribution in the reference sample having the 
highest value of feature measure |i<'Ml|. 




Figure 14. A distribution in the reference sample having a fea- 
ture that, according to the |FM1| measure, would have been as 
noticeable as the feature shown in the data 
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Class A histogram ranking among samples of 1000 histograms 



Model 


FMl 


\FM1\ 


FM2 


FM3 


Shuffled sample 


7.5 


16 


1 


1 


Simulated sample 


6.5 


10 


10 


1 



Table 1. Rankings of Feature Measures of the Class A filament 
histograms sample relative to histograms from two randomised 
samples: a sample in which the actual data is shuffled and a sam- 
ple that is constructed from a simple model. 



of the model selection function leading to figure [TT] where 
we have chose a sharp cut-off rather than a gradual one. 

There is good reason for using two different comparison 
samples. In the first sample, in which the Class A data is 
shuffled, there may be data selection issues that might them- 
selves give rise to an anomalous distribution. On the other 
hand, the artificially constructed sample, while random by 
construction, might not accurately reflect such anomalies. 
The agreement between the two samples is encouraging and 
argues strongly that there is no anomaly arising from, say, 
the filament finding process. 

It will also be noticed that in both figures it looks like 
there may be an excess of galaxies having spin angles with 
cos 6 > 0.8. That is one of several alternative explanations 
for the observed distributions. We shall quantify this below 
using a feature measure designed to explore this feature. 

T.2 Results 

7.2.1 FMl and \FM1\ 

Larger values of the feature measure \FM2\ indicate a sig- 
nificant hump or hole in the distribution of cos 6 relative 
to our reference sample. The distribution of this statistic is 
shown in figure [T^ which shows that the chances of getting 
a feature as remarkable as the one in our sample is on the or- 
der of 1%. The mean value of |_FM1| is 6.37, whlie the value 
for our Class A sample is 21. In particular, from FM2, the 
chance of getting a feature having the same shape as the 
data suggests is < 0.5%. 

We show two examples of high values of |FM1|. In fig- 
ure [13] we show the distribution having the highest value of 
FMl = 23. The Class A sample has FMlciassA = 21. The 
other example, figur JT^ shows an equally impressive feature 
having the opposite sign for FAIL Such a distribution would 
have been considered remarkable and has to be taken into 
account when assessing how significant our result might be. 

7.2.2 FM2 

The feature measure FM2 can reveal an excess of high cos 6 
objects relative to the reference model. The SDSS sample 
of Class A filaments is ranked joint 8th on FM2 along with 
6 others. The expected value for FM2 is {FM2) = 29 for 
a sample of this size, while the value for our Class A fila- 
ments is FM2ciassA = 39. The high cose part of the SDSS 
sample appears to be unusually high (probability < 1%). 
Of course that might be expected since if there is a hollow 
in the distribution relative to the model, there has to be a 
compensating excess somewhere else. 
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Figure 15. Two examples from the reference distributions having 
equal or higher values of the feature measure FM2 than the Class 
A sample. One sample was selected because it has the highest 
value of FMl in that set, the other because it has the highest 
FM2 with a positive FMl 
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Figure 16. The best of the reference models as gauged from the 
FM3 feature measure. The similarity is remarkable: this attests 
to the effectiveness of using feature measures to assess this kind 
of data. 



We show in figure [15] examples of distributions of dis- 
tributions having equal or higher i^Af2-ranked values. The 
fiatness of the region where cos 9 > 0.7 in these simulations 
is not surprising since the distribution would be flat for a 
uniform distribution of orientations. Censoring the angle rj 
that the filament makes with the line of sight removes fil- 
aments in such a way as to approximately preserve the ex- 
pected flatness of the cos 6 distribution for larger values of 
cos 9. 

We note from the sample shown in figure [TJ] having the 
highest value of of FMl, that the low cos 9 distribution is 
weaker than in our Class A sample. It is because of that 
that the composite index FA'13 for our Class A sample is 
considerably higher than for any sample in the reference 
distribution. 
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7.2.3 FM3 

FA'13 is a composite measure assessing the entire distribu- 
tion from the point of view of features. The Class A SDSS 
sample ranks first on this measure by a long way: this is 
perhaps not surprising since the Class A sample is in the 
top 1% of both FMl and FM2. However, FM3 aUows us 
to select the "best" simulation from the reference sample, 
this is shown in figure 1161 The top FM3 samples include 
a variety of values of FMl and FM2, so the selection of 
this particular example was taken from the highest ranking 
FM3 objects in such a way as to rank highly in both FMl 
and FM2. This is a rare example among the reference dis- 
tributions and if, on that basis, we were to assess the chance 
of getting this at random it would be on the order of 0.1%. 
However, any single instance of a general distribution is, by 
its very nature, rare, and we should not take that too seri- 
ously! 



7.3 Other statistical measures 

The Feature measures are somewhat ad hoc, and that is why 
they are a measure on "interestingness" rather than a mea- 
sure of "significance". We used a data mining approach to 
constantly pare down the sample to something useful with- 
out actually selecting for the very effect we were seeking to 
establish. 



7.3.1 Significance tests 

The white histogram of figure [10] shows the result of as- 
signing the sample galaxies to a randomised sample of fil- 
aments. The x-squared difference between this randomised 
histogram and the data (red) is 26.5 for 8 degrees of free- 
dom. The actual sample is very unlikely (P < 0.1%) to have 
been drawn from a sample of galaxies that are randomly 
oriented relative to the parent filament. There are many re- 
arrangements of the histogram that would lead to as bad a 
chi-squared value: what the feature measure attempts to do 
is to single out the interesting rearrangements. 

This situation arises because the x-square test takes 
no account of the distribution of the deviations: the under- 
lying hypothesis is that the deviations are independently 
distributed random errors. This is beautifully illustrated in 
Anscombe's incisive discussion in which he presents four 
quite different data samples distributions, Ansco mbe's Quar- 
tet, h aving the same basic statistical properties (jAnscombel 
I1973D . The Feature Measure is simply a way of carrying out 
Anscombe's admonishment for exploratory data analysis. 

An alternative is to ask what is the probability that the 
Class A sample as shown in figure [6] is drawn at random 
from the host sample of all 492 edge-on galaxies of figure 
[8] Applying the Mann- Whitney rank sum test to the distri- 
bution of cos 9 in the Class A sample and the total sample 
shows that the Class A sample is almost certainly not drawn 
randomly from the total sample (confidence level of the U- 
statistic — 6.9a). 



7.4 Further remarks 

7.4.1 Sample selection 

Most of the galaxies in the sample do not have b/a < 0.2 
and that automatically gets rid of a large fraction of galaxies 
in the DR5 catalogue, leaving some 50,000. However, most 
of these 50,000 or so edge-on galaxies do not lie on MMF 
filaments. That reduces the number for this analysis to below 
500: this is the sample of edge-on galaxies that lie in MMF 
generated filaments. We could have acquired a larger same 
by using a less restrictive axial ratio cut-off, like b/a < 0.25, 
but this might have introduced further uncertainties. 

The classification that was done visually afterwards 
merely defined the subset of filaments that were considered, 
via a variety of criteria, to be reliable. It turned out, with 
subsequent analysis, that this final selection favoured fila- 
ments that were both oriented some way from the line of 
sight and that were not at the extreme depths of the survey. 

The question might arise as to the accuracy of the deter- 
minations of the various angles, especially their cosines. This 
reflects on our choice of binning. As figure [191 shows, there 
should be little doubt about the accuracy of the value of 
the orientation of the galaxy relative to the sky. The prob- 
lem is in estimating the error in assessing the inclination 
of the filament to the line of sight. There is no clear pre- 
scription for doing that, especially when the filament is not 
straight. We ameliorated problems arising from estimation 
of filament orientation by including in our filament selec- 
tion process the requirements that filaments be straight and 
that galaxies should not be located at or near to filament 
branching points. The slight relaxation of these requirements 
contributes to the difference between the Class A and Class 
B samples. 

This type of selection is in fact a fundamental part of 
data mining: get rid of the objects that probably dont pro- 
vide you with any information and look at what is likely to 
tell you something. Again, that is why the focus of our at- 
tention is "interestingness" rather than "significance" . That 
should not make the analysis any less convincing, it merely 
emphasises the notion that there is something going on here 
and we should look at a bigger and better sample. 

7.4.2 Some additional comments 

Given the data, it is clear from this analysis that the distri- 
bution of spin orientations is anomalous. All Feature mea- 
sures are highly significant (better than the 1% level). Con- 
firmation of that anomaly and establishment of the nature 
of the anomaly must await a bigger sample, but that does 
not stop us from speculating about what is going on! How- 
ever, when we come to samples like DR7, the volume of data 
will preclude visual assessment of the filaments and we will 
have to do this automatically. That is perhaps the biggest 
challenge in using DR7. 



8 THE NATURE OF THE ALIGNED OBJECTS 

What is clear is that there is a dearth of objects having spin 
orientations 0.2 < cos 8 < 0.5, and this manifests itself as 
either an excess of objects with cos 6 < 0.2 or and excess of 
objects with cos 8 > 0.5, or both. This is a relatively narrow 
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Figure 17. Colour magnitude diagram for the entire sample of 
edge-on galaxies (dots) showing the gaalxies from the Class A 
filaments (circles) and the subsample of 14 candidate galaxies 
(squares). The distributions are consistent with the hypothesis 
that the candidate sample they drawn randomly from the entire 
sample. 
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Table 2. The 14 candidate galaxies that provide our evidence 
for alignment processes in filaments. They are listed in order the 
angle the spin axis makes with the host filament. Mr and Mg 
denote the absolute magnitudes of the object in the SDSS r and 
g bands, calculated using the redshift z. The ID is the filament 
number in the F ILCAT-0 catalogue o f filaments found using the 
MMF technique llAraeon-Calvo ll2007h . All of these objects, with 
the exception of galaxy 101, look like normal edge-on galaxies. 
However the following should be noted: "tidally disturbed, ''the 
SDSS published diameter looks wrong so this datum is uncertain, 
'^has a close-by neighbour of unknown redshift, "^not an edge-on 
galaxy. 
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Figure 18. The diameter - magnitude relationship for two sub- 
samples of edge-on galaxies taken from Class A filaments. The red 
dots are the sample of 14 objects having spin axis - filament angle 
such that cos 6 < 0.2, the blue boxes are a sample of galaxies 
having cos6 > 0.8. 



range of angles (60° < 6 < 80°), though it is difficult to 
make this range more precise with the data as it stands. 

The 14 anomalous objects are tabulated in table [2] and 
an instance from this table is shown in figure UM 

The colour-magnitude plot for the entire sample of edge- 
on galaxies in the sample and for the candidate sample is 
shown in figure 1171 There is no evidence on the basis of 
colours or spectra that these galaxies are in any way special 
or peculiar. It should be noted that we have not made any 
attempt to separate out early and late type galaxies. 

The diameter - magnitude relationship is shown in fig- 



ure [18] for two sub-samples of galaxies: one whose spin axes 
are perpendicular to their host filament, cos 8 < 0.2, and 
another whose spin axes lie along the filament, cos 9 > 0.8 
. Although the evidence is not strong, there is a suggestion 
that there is a lack of the bigger and the intrinsically brighter 
objects among the 14 "specials" we have identified. If this 
were true it would be an important clue in coming to an 
understanding of this phenomenon. 



9 CONCLUSIONS AND DISCUSSION 

We have taken a substantially different approach to the 
problem of finding evidence for non-randomness in the ori- 
entation of galaxies. The process we have used has been al- 
most archaeological or forensic in nature in that we look for 
likely spots for evidence of tidal influence and then search 
for objects of a specific kind. 

Firstly we have chosen to look at the relationship be- 
tween the spin axis of the galaxies and the filamentary struc- 
tures that host them. In so doing we are looking where we 
expect to find evidence. A clustered environment seems less 
promising since galaxy interactions and orbital dynamics 
will have played a role, whereas we would not expect to see 
appreciable local dynamical influences in filaments. More- 
over, the filaments are relatively easy to find using tech- 
niques like MMF based on DTFE density sampling. Indeed, 
we would argue that using a less effective filament finder 
might hinder the discovery any such alignment effects. 

The filament data was pre-processed so as to correct 
for the finger-of-god effects, and then thinned so as to bring 
the member galaxies closer to the spine of the filament. 
Then the catalogue of filaments was censored so as to give 
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a smaller, better defined sample of filaments, but without 
regard to the orientation of the galaxies lying in the fila- 
ments. In making this selection we studied the finger of god 
corrected maps in which the distribution had been thinned, 
and subjectively eliminated objects that were complex, hav- 
ing branching structures. This reduced the original sample 
of 426 filaments to 67 "Class A" filaments. Most contained 
but one edge-on galaxy: only three pairs of our sample of 69 
galaxies in Class A filaments lie in the same filament and no 
filament contains more than one of our "special" subgroup 
of 14 galaxies. No Class A filament contained more than two 
edge on galaxies. 

We suggest that the distribution is not consistent with 
a random distribution of spin orientations and, in particular, 
we find 14 objects among these 67 filaments that are oriented 
perpendicular to their host filament. There is corroborating 
evidence for this assertion, The slightly wider sample of 121 
Class A and Class B objects gave some support for this, as 
does the sub-sample of filaments lying in the plane of the 
sky. 

In order to assess the statistical standing of this con- 
clusion we characterise the distribution of spin angles by 
a number of "feature measures" and assess the chances of 
coming up with values for those measures equal to or greater 
than the values for our Class A SDSS sample. Both shuffling 



the spins in the Class A sample among their filaments and 
using a sample of 1000 histograms for randomly distributed 
orientations shows that our observed Class A distribution is 
indeed remarkable. 

If we take this analysis at face value, there are at least 
two possible interpretations of the histogram we find for 
the spin orientations. Either there is an excess of objects 
having their spin axes perpendicular to the spine of the host 
filament, or there is an excess of filaments having their spin 
axes parallel to the spine. These views are not inconsistent: 
it could simply be that the distribution of orientations is 
bimodal. It is easier to follow up the first of these alternatives 
since there are only 14 objects involved, while in the second 
of the alternatives there is no way of deciding which galaxies 
might individually be responsible for the excess. 

The distribution of the relative orientations of the spin 
axis of the galaxy and its host filament provided us with a 
sample of 14 candidate galaxies that we would claim pro- 
vide evidence for such tidal interactions. These candidates 
are typical of any other edge-on galaxy in the entire sample 
- they are not in any way distinguished except by their un- 
expected orientation and perhaps in being slightly smaller 
than other edge-on galaxies in the sample. 

It turns out that these candidate galaxies lie centrally 
in relatively insignificant filaments. This may suggest that 
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they have suffered less from interactions with other galaxies 
and have been influenced mostly by the global tidal field 
exerted by their parent filament. 

ft might occasion surprise that an original sample of 
over 50,000 galaxies has been whittled down to a mere 14! 
The first sample of filaments was defined using the entire 
galaxy sample, and we selected only those containing edge- 
on galaxies defined as having axial ratio b/a < 0.2. This 
left us with only 426 filaments, which subsequent quality 
assessment reduced to 67 containing 69 edge-on galaxies. 

Obviously this should now be repeated using a sam- 
ple like the DR7 release of the SDSS. If these 14 galaxies 
are indeed fossils of earlier tidal interactions, using a bigger 
catalogue would either give many more candidate objects 
or eliminate the effect entirely. One advantage of a bigger 
sample is that we could make a fiatter axial ratio cut. The 
general procedure is, however, far less straightforward since 
it would not be possible to do the visual triage of filaments 
to select out "clean". Class A, filaments. This would have 
to be done algorithmically. 

However, until that is done, we offer this study as a set 
of fossils providing forensic evidence that tidal interactions 
have affected the orientations of the spin axes of galaxies. 
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APPENDIX A: THE TIDAL SHEAR FIELD 

The relation between the fluctuating component of the den- 
sity field, 5{r,t), and gravitational potential "I> is established 
through the Poisson-Newton equation: 

= 47rGp„(t)a(t)2 S{r,t), (Al) 
5(r,t) = ^(£l^. (A2) 

pu 

Here pmit) is the mean density of the mass in the universe 
that can cluster (dark matter and baryons). 

The peculiar gravitational acceleration due to the in- 
tegrated effect of all matter fluctuations in the Universe is 
related to $(r, t) through g — — V$/a: 

g(r, t) = -4^Gp„(t)a(t) / dr' 5{v' , t) i^— . (A3) 

and drives the peculiar (non-Hubble expansion components) 
of the cosmic motion. The cosmological density parameter 
0,m{t) is defined by pu, via the relation QmH^ = (87rG/3)pm 
in terms of the Hubble parameter H. 

The tidal shear arising from this accelera- 
tion is given by the (traceless) tidal tensor Tij 
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l|van de Weygaert &: Bertschingerlll996l ), 



1 J dgb,i ^ dgb.j 



2a 1 dxi dxj 



+ ^(V.g.)5, 



(A4) 



This can be expressed directly in terms of the fluctuation 
component of the density field, 5, via the equation 

T^J{r) = j 5(r')Q.,dr - \nH^ S{r,t) (A5) 

3(ri — ri)(r'j — rj) — \r' — S 



(A6) 



This expression shows explicitly that source of the tidal 
field is the quadrupole component of the fiuctuating mat- 
ter distribution, QijS{r'). The quadrupole component of 
the field falls of as r~^, so, unlike the isotropic component, 
it is localised. However, there is a considerable contribu- 
tion from the larg est scales as can be seen in the study of 
void alignments bv lPlaten et al. I (l2008l ) and in the study of 
shear fields in relatio n to gravitational lensing (|Gunn II 19671 : 
IWittam et"ani2000l '). This can produce systemic shear that 
is correlated over large scales, resulting in systemic align- 
ments that survive until mergers and tidal forces from other 
nearby structures start playing a role. 

The tidal shear tensor has been the source of intense 
study by the gravitational lensing community since it is 
now possible to map the distribution of large scale cosmic 
shear using weak lensing d ata. Examples are the studies by 
iHirata fc SeTia^ (|2004l ) and lMassev et al. I (|2007l ). 



APPENDIX B: TIDAL ALIGNMENT OF 
GALAXY SPINS 

One may obtain insight into the issue of the alignment of 
a halo's angular momentum L with the surrounding mat- 
ter distribution by evaluating the expression for the angular 
momentum L, 



(Bl) 



(with eijk the Levi-Civita symbol) in the principal axis frame 
of the tidal tensor Ta, 
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(B2) 
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— A2) I12 , 





where (Ai,A2,A3) are the tidal field eigenvalues. Since, by 
definition, (A3 — Ai) is the largest coefficient, in a statistical 
sample of halos L2 will on average get the largest contribu- 
tion if we assume that the tidal and inertial tensors Tij and 
lij are uncorrelated. 

However, the assumption above is not well just ified, and 
in fac t the two t ensors are f ound to be correlated l|Porcianil 
|2002| ). Recently, iLee et al] l|2009l ) numerically determined 
the two-point correlations of the three eigenvalues of the 
nonlinear traceless tidal field in the frame of the principal 
axes of the tidal field. The numerical findings indicate that 
the correlation functions of the traceless tidal field and the 
density field are all anisotropic relative to the principal axes. 
Interesting is also that their correlations have much larger 
correlation length scales than that of the density field and 



increase along the directions normal to the first principal 
axes of the tidal field. 

To incorporate the effect of a c orrelation b etween the 
tidal and inertial tensors Tij and Jfci. lLee fc PenI (2000 ) and 
I Lee fc PenI (|200ll ) suggested an analytical formulation in- 
volving a useful parametrization of the correlation. This they 
accomplished by writing down an equation for the autocor- 
relation tensor of the angular momentum vector in a given 
tidal field, averaging over all orientations and magnitudes 
of the inertia tensor. Arguing that the isotropy of the un- 
derlying density distribution allows the replacement of the 
statistical quantity {Imqins) by a sum of Kronecker deltas, 
one is left with the result that the autocorrelation tensor of 
the angular momentum vector is given by 

{UL,\T) OC + (i^,, - T.kTkj) (B3) 

If it is asserted that the moment of inertia and tidal shear 
tensors were uncorrelated, we would have only the first term 
on the right hand side, ^Sij-. the angular momentum vec- 
tor would be isotropically distributed relative to the tidal 
tensor and one would not expect the presence of any signif- 
icant galaxy spin alignments. Recognizing that the inertia 
and tidal t e nsors in gen eral are not mutually independent, 
iLee fc PenI (|20o3,[20o3) introduced a parameter c to quan- 
tify this correlation. 



OC -Oi 



kj) 



(B4) 



where c = for randomly distributed angular momentum 
vectors. The case of mutually dependent tidal and inertia 
tensors is described by c = 1 (see equation IB3|I . Finally, 
they introduce a different parameter a = 3c/5 and write 

^ ^ -<5ij - aTikTkj (B5) 



{LiLj\T) OC 



which forms the basis of much current research in this field. 
The value der ived from a recen t study of the Millennium 
simulations bvlLee fc PenI (1200711 is a « 0.1. 

iLee fc Erdogdu I (2007) quantified the preferential 
alignment of the angular momentum vector L along the in- 
termediate 

principal axis T2 of the tidal tensor, by deriving the 
corresponding probability distribution p{92) for the align- 
ment angle cos 62 = |L ■ T2|/|L|. Including the effect of the 
correlation between the tidal and inertial tensor by means 
of the parameter c, the find that 



p(cos t 



(1 + c) 



l-£ 
2 



1 



c 11 — - cos 6*2 



3/2 ' 



(B6) 



where cos 62 is assumed to be in the range [0, 1]. If c = 0, 
the probability distribution will be a uniform distribution 
pi.02) = 1. 



APPENDIX C: FINDING FILAMENTS WITH 
MMF 

One of the key issues in this investigation is to rigorously 
identify filaments in the large scale structure. The Multiscale 
Morphology Filter, MMF, that we use to identify filamen- 
tary structures is presented in detail in (|Arag6n-Calvo et al.l 
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I2007I ) where a detailed step-by-step description of the MMF 
algorithm can be found. In this Appendix we briefly sum- 
marize the steps involved in the morphological segmentation 
of the cosmic web obtained from the N-body cosmological 
simulation. 



Structure 
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A constraints 
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Al 
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CI Scale Space 

The DTFE density field /dtfe is the starting point of the 
morphological segmentation. The density field is smoothed 
over a range of scales by means of a hierarchy of spherically 
symmetric Gaussian filters Wg having different widths 
The n*'' level smoothed version of the DTFE reconstructed 
field /dtfe is assigned /„, 



where Wg denotes a Gaussian filter of width _R„: 



WG{y,x) = 



(27ri?2)3/2 



exp 



\y - x[ 

2Rn 



(CI) 



(C2) 



Scale Space itself is constructed by stacking these vari- 
ously smoothed data sets, yielding the family <1> of smoothed 
density maps /„: 



$ = 



U 



(C3) 



Table CI. Behavior of the eigenvalues for the characteristic mor- 
phologies. The lambda conditions describe objects with intensity 
higher that their background (as clusters, filaments and walls). 
From the constraints imposed by the A conditions we can de- 
scribe the blob morphology as a subset of the line which is itself 
a subset of the wall. 



\a{x) and arrange them so that Ai > A2 > A3: 



with 



0, a = 1,2, 3 



Al > A2 > A3 



(C5) 



The \i(x) are coordinate independent descriptors of the be- 
haviour of the density field in the locality of the point x and 
can be combined to create a variety of morphological indi- 
cators. The criteria we used for identifying a local bloblike, 
filamentary or sheetlike morphology are listed in table [C2l 



A data point can be viewed at any of the scales where scaled 
data has been generated. The crux of the concept of Scale 
Space is that the neighbourhood of a given point will look 
different at each scale. There are potentially many ways of 
making a comparison of the scale dependence of local envi- 
ronment. We address the local "shape" of the density field. 



C2 Local Shape 

The local shape of the density field at any of the scales _R„ 
in the Scale Space representation of the density field can 
be quantified on the basis of the Hessian matrix, T-lij — 



V«j/n(x), 



dxidx 



■Mx) = /di 



dxidx 



■Wg{Ru 



{xj - yi){xj - yj) - 5^jRl 



WGiy,x) (C4) 



where {xi, a:;2, 2:3} = {x, y, z} and Sij is the Kronecker delta. 
In other words, at each level n of the scale space representa- 
tion the Hessian matrix is evaluated by means of a convolu- 
tion with the second derivatives of the Gaussian filter, also 
known as the Marr (or, less appropriately, "Mexican Hat") 
Wavelet. In order to properly compare the values of the Hes- 
sian arising from the differently scaled variants of the data 
that make up the Scale Space, the Hessian is renormalized, 
a — R^y., where Rs is the filter width that has been used. 

The eigenvalues Xi of the Hessian matrix determine the 
local morphological signal, dictated by the local shape of 
the density distribution. A small eigenvalue indicates a low 
rate of change of the field values in the corresponding eigen- 
direction, and vice versa. We denote these eigenvalues by 



C3 Multiscale Structure Identification 

In practice, we are interested in the local morphology as a 
function of scale. In order to establish how it changes with 
scale, we evaluate the eigenvalues and eigenvectors of the 
renormalised Hessian "H of each dataset in the Scale Space 

Since we are looking for three distinct structural mor- 
phologies - blobs, walls and filaments - the practical imple- 
mentation of the segmentation consists of a sequence of three 
stages. Because curvature components are used as structural 
indicators, the blobs need to be eliminated before looking for 
filaments, after which the filaments have to be eliminated 
before looking for walls. This results in the MMF procedure 
following the sequence "clusters — >■ filaments — >■ walls". At 
each of these three steps, the regions and scales are identi- 
fied at which the local matter distribution follows the cor- 
responding eigenvalue signature. 

In practice, the MMF defines a set of morphology 
masks, morphology response filters and morphology filters 
for each of the three different morphological components: 
clusters, filaments and walls. Their form is dictated by the 
particular morphological feature they seek to extract, via 
the eigenvalues at each level in scale space and the crite- 
ria for each of the corresponding morphologies (table [C2)) . 
The local value of the the morphology response depends on 
the local shape and spatial coherence of the density field. 
The morphology signal at each location is then defined to 
be the one with the maximum response across the full range 
of smoothing scales. 

The end result is a map segmented in clusters, filaments 
and walls that have been identified as most outstanding, and 
which vary in scale over the full range of scales represented 
in the Scale Space. 
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Figure Dl. Distribution of spins in filaments with censored distribution of angles to line of sight. See also figure[4] The points A, B, C, D 
lie on the base of a circular cone of semi vertex angle 9. If there were no censorship of the filament orientation, all angles 9 would be 
equally accessible. However, if we eliminate filaments making angles less than rjc with the line of sight, the sample of spin orientations 
9 > 7r/2 — rjc will be truncated. Only spin axes at angle 9 to the filament lying along the segments AB and CD will be in the sample. 
Since the distribution of 9 values around ABCD would be uniform the fraction of objects in the censored sample is the ration of the 
lengths of the segments AB and CD to the circumference of the circle. 



APPENDIX D: DISTRIBUTION OF RANDOM 
SPINS WITH CENSORED FILAMENT 
ORIENTATIONS 

The angle 6 between the spin vector of a galaxy and its host 
filament is 



I cos Q sin < 



(Dl) 



(see figure 3] for the definitions of the angles and equation 
Q). If the filament orientation is random and spins direc- 
tions are uncorrelated with the host filament, then 6 is ran- 
domly distributed over the unit sphere and cos 6 is uniformly 
distributed on the interval [0, 1] . 

We have seen that, in our catalogue of filaments, the 
high-grade Class B and Class B filaments tend to lie trans- 
verse to the line of sight (see figure [5]). This imposes a con- 
straint on the distribution of angles 6 that the spin axis of 
a galaxy in our catalogue can make with its host filament. 
Even if the filaments and spins were randomly oriented, we 
would see systematically fewer angles 6 having low values of 
cos 6. The situation is illustrated in figure iDl] 



Dl Definitions 

The geometry of figure [DTI is as follows. The circle on which 
the points A, B, C, D lie is the base of a right circular cone 
with apex at O. The side of the cone can be taken to be of 



unit length. The axis of the cone, OR, is the spin vector of 
a galaxy in the catalogue. The point R lies in the plane of 
the base of the cone, on the axis of the cone. Thus a point P 
on this circumference represents a sample filament, and the 
vector OP makes an angle 77 with the line of sight, and an 
angle 6 with the spin axis of the galaxy. 6 is the semi vertex 
angle of the cone. 

The circumference of the base of the cone represents all 
filaments that are at angle 9 with the spin axis OR. Let us 
focus attention on the locus of points P: the circumference 
of the base of this cone. The position of P on the circum- 
ference can be parametrised by an angle a in the plane of 
the base of the cone. The line PQ in the plane of the base 
of the cone is parallel to the line of sight, and the angle a is 
measured from the line RQ (in a clockwise sense) . If every- 
thing were randomly and uniformly distributed on a sphere, 
the distribution of the angle a would be uniform. 



D2 Censorship 

The angle a is seen to be given in terms of 77 and 9 by the 
simple relationship 



cos r] 
sin 9 



(D2) 
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Figure D2. Comparison of Monte Carlo simulation of randomly 
oriented galaxies selected from a censored sample of randomly 
oriented filaments. The censorship angle is eta > 45° , as in figure 
111! The Monte Carlo simulation is shown as the light grey his- 
togram. The dark grey histogram is equation IIDSI l appropriately 
normalised. 



This follows because PQ = cost; and PR — sin 6*. From this 
we have that the values of a traced out as ri varies for a 
given 9 are given by 



a — arcsm 



cos Tj 



cos?7 < sinO. (D3) 



The constraint cos rj < sin 6 shows that not all 0- values are 
attained at all values of rj. For a given value of 9 all values 
of a are attained only if 

r] > Tjc where cos77c = sin^. (D4) 



Geometrically, this is depicted in figure [DTI where the base 
of the cone is shown as the (projected) circle containing the 
points A, B, C, D and P. 

Consider now a constraint such that our catalogue con- 
tains only filaments with rj > rjc, for some cutoff rjc- The 
angle rj, as depicted in the figure, is large enough that only 
points on the segments AB and CD can be reached with 
values of rj greater than the censorship angle rjc- Points on 
the segments AD and BC have rj < etUc, and so the cor- 
responding values of a are unavailable: censorship in eta 
imposes censorship on a. 



D3 Statistical distribution for random 
orientations 

In the absence of any r;-censorship the distribution would be 
uniform in cosO. When a cut-off, rjc, in rj is introduced, the 
uniform distribution is attenuated by a factor 



2 2 
/ — —a = — arcsin 

TT TT 



COS rjc 



6 > 11/2- jjc 
otherwise. (D5) 



This can easily be tested using simulations, one of which 
has been displayed in the grey histogram of figure 1111 The 
comparison between equation (jDSjl and figure [TT] is shown 
in figure lD2l 
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